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Abstract: The first stars to form in the Universe may be powered by the annihilation 
of weakly interacting dark matter particles. These so-called dark stars, if observed, may 
give us a clue about the nature of dark matter. Here we examine which models for particle 
dark matter satisfy the conditions for the formation of dark stars. We find that in general 
models with thermal dark matter lead to the formation of dark stars, with few notable 
exceptions: heavy neutralinos in the presence of coannihilations, annihilations that are 
resonant at dark matter freeze-out but not in dark stars, some models of neutrinophilic 
dark matter annihilating into neutrinos only and lighter than about 50 GeV. In particular, 
we find that a thermal DM candidate in standard Cosmology always forms a dark star as 
long as its mass is heavier than ~ 50 GeV and the thermal average of its annihilation cross 
section is the same at the decoupling temperature and during the dark star formation, 
as for instance in the case of an annihilation cross section with a non-vanishing s-wave 
contribution. 
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1. Introduction 

The first stars, also referred to as Population III stars, are the first luminous objects in the 
Universe. They contribute to the reionization of the interstellar medium, they provide the 
heavy elements (metals) that eventually become part of the later generations of stars, and 
they may be the seeds of the very massive black holes observed in quasars. 

It was shown in III |2], [| that the first stars to form in the Universe may be powered by 
the annihilation of dark matter particles instead of nuclear fusion. These dark-matter pow- 
ered stars, or dark stars for short, constitute a new phase of stellar evolution. Besides the 
assumption that dark matter is made of weakly interacting massive particles (WIMPs) that 
can self- annihilate into ordinary particles, three conditions are necessary for the formation 
of a dark star. 

The first condition is that the density of dark matter at the location of the (proto)star 
must be high enough for dark matter to efficiently and rapidly annihilate into ordinary 
particles, releasing a large amount of energy. The first stars are believed to form at the 
center of dark matter halos when the Universe was young (redshift z ~ 10-50) and denser 
than today. Not only the dark matter density at the center of those early halos was high, 
but as the baryonic gas contracted into the first protostars, more dark matter was gathered 
around the forming object by the deepening of the gravitational potential (gravitational 
contraction). Cosmological parameters and the evolution of the gas density completely 
determine the resulting density of dark matter at the location of the first protostars. An- 
alytic and numerical evaluations [2|, |]] lead to a resulting density which is high enough 
to satisfy the first condition for the formation of a dark star. 



The second condition is that a large fraction of the energy released in the dark matter 
annihilation must be absorbed in the gas that constitutes the (proto)star. The fraction 
/q of annihilation energy deposited into the gas depends on the nature of the annihilation 
products. Typical products of WIMP annihilation are charged leptons, neutrinos, hadrons, 
photons, W and/or Z bosons, and Higgs bosons. The latter (W, Z, and Higgs) decay 
rapidly into leptons and hadrons. The hadrons themselves, which are mostly charged and 
neutral pions) decay rapidly into charged leptons, neutrinos, and photons (although a small 
number of stable particles like protons can also be produced). After ~ 10 -8 seconds, all 
unstable elementary particles, including the muon, have decayed away, and only protons, 
electrons, photons and neutrinos survive. Protons have a large scattering cross section 
with the protostar medium and are quickly absorbed. Electrons and photons can ionize 
the medium and/or generate electromagnetic showers. For WIMPs with mass m > 0.5 
GeV, electromagnetic showers are the dominant process. At the time when the following 
third condition for a dark star is satisfied, the protostar has a diameter of more than 40 
radiation lengths, implying that all the energy released in protons, electrons, and photons 
is absorbed inside the protostar. Only the fraction of energy carried away by the neutrinos 
is lost for what concerns a dark star. 

The third condition for the formation of a dark star is that the heating of the (proto)star 
gas arising from the dark matter annihilation energy must dominate over any cooling 
mechanism that affects the evolution of the (proto)star. In ||], it was shown that the dark 
matter heating rate Qdm> m energy deposited per unit time and unit volume, is given by 
the expression 

n _ , Hds 2 /, U 

Ql)M - JQ P , (1-1) 

m 

where p is the dark matter density inside the (proto)star, which is determined by the 
cosmological model, and (crv)ds is the average value of the dark matter annihilation cross 
section a times WIMP relative velocity v inside a dark star. To the extent that electro- 
magnetic showers are generated, i.e. m > 0.5 GeV, all dark star properties depend on the 
particle physics model only through the quantity 

/Q^- (1-2) 
m 

Ref. Q fixed the annihilation cross section to (&v)ds = 3 x 10~ 26 cm 3 /s and examined 
a range of WIMP masses m from 1 GeV to 10 TeV. In addition, Ref. [Q] assumed Jq = 2/3, 
based on simulations of neutralino dark matter annihilation in the Minimal Super symmetric 
Standard Model (MSSM). For this range of Qdm, they compared the heating and cooling 
rates along protostar evolution tracks from [0], and concluded that there is a time during 
the evolution of the protostar in which the dark matter heating dominates over all cooling 
rates. This finding lead to the realization that dark stars may be possible. 

In this paper, we examine the possible values of Qdm for a large selection of particle 
physics models, and verify if the third condition above is satisfied in these models. We 
find that not all particle dark matter models lead to the formation of dark stars, although 
the models that do not form dark stars are either tuned to resonant annihilation or rather 
artificial. 
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Figure 1: Condition for the formation of a dark star, in terms of the protostar gas density and 
temperature. The gray band shows possible evolution tracks of the protostar obtained through 
numerical simulations in a ACDM cosmology ||. The red lines show critical curves on which 
the heating rate from dark matter annihilation equals the total cooling rate of the protostar gas. 
Critical curves are labeled by the value of fQ(&v) c i s /m in units of cm 3 s _1 GeV -1 . A dark star 
forms at the intersection of a critical line with the gas evolution track. No dark star can form for 
fQ(<Jv)ds/m < 10 -32 cm 3 s _1 GeV -1 (solid critical line on the right). 
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The restriction imposed by the third condition for dark star formation is best expressed 
in terms of a condition on the quantity fQ{crv)d s / m - Following we have computed 
the critical lines in the gas temperature-density plane at which the heating rate from 
dark matter annihilation equals the total cooling rate. These lines are shown in Figure [j] 
for a wide range of values of fQ(<Jv)ds/ m , from 10~ 18 cm 3 s" 1 GeV" 1 to 10~ 32 cm 3 s" 1 
GeV -1 . Below the latter value, the heating-cooling critical line no longer intersects the 
thermodynamic track of the protostellar gas, indicated by the gray band obtained through 
numerical simulations of the formation of the first stars in a ACDM cosmology || . In other 
words, for fQ(uv)ds/ m < 10~ 32 cm 3 s" 1 GeV , the protostar is expected to contract to 
a regular Population III star powered by nuclear fusion without passing through the dark 
star phase. At the other side of the fQ{crv)ds/ m range, the critical line reaches a limiting 
curve given by the vertical line labeled fc){crv) ds/m = 10~ 18 cm 3 s" 1 GeV -1 . Larger values 
of fQ((Jv)ds/m give the same vertical line. Thus, as expected, if the annihilation rate is 
large, a protostar passes through the dark star phase. Therefore the third condition for 
the formation of a dark star is 



(av) 



ds 



> 1 x 10~ 32 cm 3 s -1 GeV" 



m 

3 x 10~ 26 — 3 
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The choice (cf)ds = 3 x 10" 2e cm 3 /s in |lj was motivated by the assumption that 
the dark matter WIMPs are produced thermally in the early Universe. That is, that 
the WIMPs are generated in matter-antimatter collisions at temperatures higher than 
Tf ~ to/20, which is the temperature after which WIMP production "freezes out" and 
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the comoving WIMP number density remains (approximately) constant. Ref. [|| used the 
following simple inverse-proportionality relation between the present WIMP density £l x 
and the annihilation cross section (av)f Q at the time of WIMP freeze-out, 

n , 2 3 x l(T 27 cm 3 /s 

n x h = 7 — ; — ■ L4 ) 

Furthermore, ref. ||] simply assumed that the velocity-averaged annihilation cross section 
times relative velocities at the time of freeze-out and in a dark star have the same value, 
(o"f)ds = ( " u )fo- In reality, the relation between Q x and {av)f Q is more complex, and in 
addition (av )d s may differ from (av)f because av may depend sensitively on the WIMP 
velocity v. In this regard, we notice that the average WIMP speed at freeze-out is of the 
order of 

Ufa ~ \ — f= ~ 7 x 10 4 km/s, (1.5) 

V m ^20 

while the typical speed of WIMPs in a dark star can be estimated from the orbital velocity 



GM , . , , 

Vds ~ y — 30 to 300 km/s, (1.6) 

namely ~ 30 km/s for a newly-born 1-Mq dark star of 1 AU radius or ~ 300 km/s for a 
mature 600- M & dark star of 5 AU radius. 

A neutralino in the MSSM provides an example of a more complex relation between 
U x and (cru)fo. At the same time, it allows the direct evaluation of both (av)d s and /q, 
and in general it has (av)d s / (crv){ . Section |2| explores this case. 

Kaluza-Klein dark matter is examined in Section 3, where it is concluded that generi- 
cally in these models (av)ds tends to be larger or comparable to {av)f . 

Leptophilic models of dark matter proposed to explain the PAMELA positron excess 
and the Fermi and HESS cosmic-ray electron-positron data provide another example in 
which (av) <j s may not be the same as (crv)f Q . They are examined in Section ||. 

Finally, we push /q down using dark matter particles that annihilate exclusively into 
neutrinos ( "neutrinophilic" models). In these models, even if annihilation produces pre- 
dominantly neutrinos that escape the forming star, W- and Z-bremsstrahlung processes 
may generate enough charged leptons to actually form a dark star. We examine this case 
in Section || 



2. MSSM 

Because of its many free parameters (more than 100), the MSSM provides a variety of ex- 
amples in which the annihilation cross section in the dark star differs from the annihilation 
cross section at the time of freeze-out, or more precisely (av)^ ^ {crv)f . 

There are several ways in which the equality {crv)d s = {^v)f can be violated in the 
MSSM ||. First, the quantity av may depend on the relative velocity v. This includes 
three cases: (i) p-wave annihilation in which av = a + bv 2 is dominated by the bv 2 term 
at freeze-out (here a and b are constants); (ii) resonant annihilation in which av follows 



a Breit-Wigner function (1 + v 2 ) c /[(v 2 + 5) 2 + 7 2 ], where 5, 7 and c are constants; and 
(iii) threshold annihilation in which an annihilation channel is kinematically accessible at 
freeze-out but not in a dark star thanks to the higher particle kinetic energies at freeze-out. 
Second, the annihilation reactions that determine the freeze-out time may be unrelated to 
the neutralino-neutralino annihilation that occurs inside a dark star, in that the freeze- 
out temperature may be high enough to convert neutralinos into heavier supersymmetric 
particles that annihilate much faster (a phenomenon called coannihilation). It is then 
the annihilation cross section of the heavier supersymmetric particles that determines the 
neutralino relic density, and this cross section is in general not the same as the neutralino- 
neutralino annihilation cross section, thus {av)ds 7^ (c^)fo- We remark in passing that 
resonant annihilation and coannihilations are not rare phenomena in the MSSM, and are 
actually essential to obtain a neutralino dark matter in the minimal supergravity or con- 
strained MSSM models. 

To illustrate these four cases (p-wave annihilation, resonant annihilation, threshold 
annihilation, and coannihilation), it is sufficient to consider a so-called effective MSSM 
(effMSSM) with eight free parameters fixed at the electroweak scale pi] ]. These parameters 
are: the CP-odd Higgs boson mass rriA, the ratio of neutral Higgs vacuum expectation 
values tan/3, the Higgs mass parameter p, the gaugino mass parameters Mi and M2, the 
slepton mass parameter m^, the squark mass parameter nig, the ratios Af/m^, A^/m,q and 
A^/rriq involving the trilinear couplings A?, Aj and A^ of the third generation of sleptons 
and squarks (the three ratios are assumed to be equal). 

We consider the parameter region of the effMSSM in which the lightest neutralino is 
the lightest supersymmetric particle and its relic density £l x is within the cosmological 
range 0.098 < U x h 2 < 0.122. In this region, we compute {av)d s as the value of av at v = 0. 
For each point in this region we also compute /q as the fraction of annihilation energy that 
does not go into neutrinos. In obtaining /q, it is safe to assume that the particle cascades 
after annihilation develop in vacuum, since muons, taus and light mesons produced in the 
annihilation decay to neutrinos before being stopped in the dark star medium. 

Figure || shows the values of the combination fQ{crv) < ± s /m x obtained in the way just 
described as a function of the neutralino mass m x . There are four classes of points: (i) the 
spread of points along the direction sloping down to the right is due to p-wave annihilation; 
(ii) the V-shaped feature at m x ~ 45 GeV is due to resonant annihilation through the Z 
boson (other resonant annihilations, through the lightest Higgs boson of mass varying from 
115 GeV to 120 GeV, are visible at m x ~ 60 GeV), (iii) the "fingers" of points dropping 
from the p-wave band of points arise from threshold coannihilation, and (iv) the shaded 
region on the right of m x ~ 100 GeV corresponds to possible coannihilation with staus (the 
dashed line shows the similar boundary for sneutrino coannihilations). These four cases 
are described in the following. 

In p-wave annihilation, the dominant contribution at freeze-out comes from the p-wave 
term bv 2 in av. The p-wave contribution to av, which is instrumental to provide the correct 
neutralino relic density, is suppressed as far as the evolution of the dark star is concerned. 
In fact, for a newly-born 1-M & dark star, one has (^ds/^fo) 2 ~ 2 x 10~ 9 . Thus in this 
case, (av)ds — a while (av)f Q ~ b(v 2 ) ~ 6/20, and in general they differ. Their exact ratio 
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Figure 2: Scatter plot of the combination fQ{cv) < is/'rn x plotted as a function of the neutralino 
mass m x for an effective eight-parameter MSSM model. Below the horizontal line a dark star cannot 
form. For resonant annihilation (the V-shaped feature at m x ~ 45 GeV and the descending points 
around m x ~ 60 GeV) and coannihilations (the shaded region and the dashed line on the right of 
m x ~ 100 GeV), the quantity /q(ctw) ds/m-x ma y be *oo smau f° r a dark star to form. For p-wave 
annihilation (the band sloping down to the right) and threshold annihilation (the various "fingers" 
of points dropping from the p-wave annihilation band), dark stars would be able to form. 



depends on the particle physics parameters contained in the coefficients a and b. In our 
effMSSM scan, p-wave annihilation gives rise to a spread in fQ{<rv)<i s /m of about one order 
of magnitude (band of points sloping down to the right in Figure ^). 

The Z resonance at m x ~ 45 GeV provides an example of resonant annihilation. The 
resonant part of the neutralino-neutralino annihilation cross section is given by 



(s 



m 



2 \2 



m 2 (s — m 2 z ) 2 + T 2 z m 2 z 



2 ■ 



(2.1) 



where /3f is the speed of the final products in units of the speed of light, and g e g contains 
the coupling constants and the mixing angles of the neutralinos and of the final particles 
involved. The velocity dependence of (<jv)z can be obtained by writing s = 4m^(l + v 2 ), 
from which one finds, neglecting the mass of the final products, 

(v 2 + 5) 2 
m x (v 2 + 5) 2 + 7 2 ' 

where 5 = 1 — m|/(4m^) and 7 = Tzmz/i^m 2 ,). On resonance, that is for 2m x = mz or 
5 = and 7 = Tz/mz = 0.0273, the velocity-averaged ((crv)z) has very different values 



(2.2) 



at freeze-out (v ~ 0.2c) and in a dark star (v ~ 0). At freeze-out, the thermal average of 
[ctv)z on resonance is, using Tf Q = m x /20, {{(Jv)z){ = 0.97 g^ s /m x . On the other hand, 
in a dark star, one has on resonance, for v = 30 km/s, {{av)z)ds = 1-3 x 10~ 13 g^ s / m?.. 
While (((7w)^)fo ~ 3 x 10~ 26 cm 3 /s to provide the correct relic density, {{ov)z)ds is thirteen 
orders of magnitude smaller. These very different values of {crv)f and (crv)^ give rise to 
the V-shaped feature in Figure Q around m x = mz /2 ~ 45 GeV. (Similar resonant features 
through the lightest Higgs boson appear superposed at m x = m#/2 ~ 60 GeV.) 

Threshold annihilation occurs when the neutralino mass is slightly smaller than half 
the total mass of the final annihilation products in a specific channel (for example, XX ~~ ^ 
WW). In this case, kinetic energy is required for the reaction to occur. This kinetic 
energy is available at the time of freeze-out thanks to the relatively high temperature of 
the neutralinos, but is not available at the lower velocities of neutralinos in a dark star. 
Therefore, the annihilation into the specific channel (xx ~~ WW in the example) occurs 
at freeze-out but not in a dark star. The cross section (av)^ s is thus smaller than (av){ . 
In Figure |2| this is illustrated by the "fingers" of points dropping from the p-wave band at 
m x ~ 80 GeV (the WW channel) and m x ~ 190 GeV (the tt channel). In neither case the 
suppression of {av)d s is severe enough to bring the points outside the parameter region in 
which dark stars can form. 

For coannihilations, the relic density is determined by an effective annihilation cross 
section {<jv) e s, which is an average of the annihilation cross sections of all reactions between 
the neutralino and the coannihilating particles. In minimal supergravity models, which are 
a subset of MSSM models, coannihilations occur in specific regions of the parameter space 
in which the stau f is very close in mass to the neutralino Xi an d with more tuning 
of the parameters when the stop t is very close in mass to x- I n the general MSSM, 
coannihilations may also occur between the lightest and second lightest neutralino, and 
between the neutralino and the chargino. 

For the sake of illustration in the context of dark stars, we focus on stau coannihilations, 
because the experimental lower bound on the stau mass (m,f > 98 GeV) is smaller than 
the lower bound on squark masses and thus the coannihilation region in parameter space 
is larger. In the case of stau coannihilations, the effective annihilation cross section is 
(approximately) 

__ (av) xx + (av) xf e -( mf - m *)/ Tfo + (av) f f e ~^ m -- m ^/ T{ ° 

\ aV < eS ~ 1 + e -(m f -m x )/T {o + e -2(m f -m x )/T io " t 2 " 3 ) 

Here (av) xx , (av) x f and (crv)ff are the total annihilation cross sections for XX ~^ anything, 
XT — > anything, and ff — > anything, respectively. Since reactions like ff -4 tt are elec- 
tromagnetic processes, (av)ff ~ a 2 /m?, which is much larger than the cross section for 
XX ~^ TT i ( av )xx ~ c?rr? T lm\. In fact, for rrif = 100 GeV (1 TeV), their ratio is approxi- 
mately (av)ff/(av) xx ~ m|/m 2 > 3 x 10 3 (3 x 10 5 ). Thus with an appropriate choice of 
the mass difference rrif — m x in Eq. ( |2.3D , one can obtain an effective annihilation cross 
section three or more orders of magnitude larger than the XX annihilation cross section, 
and a relic density three or more orders of magnitude smaller than without coannihilations. 
This argument allows us to estimate a lower limit on {crv)^ s in a dark star using just the 



annihilation cross section for \X ~> TT without having to compute the relic density in the 
presence of coannihilations. The annihilation cross section for XX ~> TT can be computed 
analytically and can be limited from below by keeping only the diagram with f exchange 
and choosing appropriate neutralino and stau mixings. In this way, we obtain 

ii \ \ ^ 7rQ;2 m T In A\ 

(°(xX^TT)v) ds > — — \. (2.4) 

62 cos* t>w fn~ 

Then we set nif = m x as appropriate for stau coannihilations. Moreover, bremsstrahlung 
(XX ~^ TT'y) gives a contribution to (av) ds that exceeds the lower limit just computed in 
Eq. (|2.4j ) at large neutralino masses. For m x = nif we estimate 

a 3 

(v(XX^TT~f)v) ds ~ —j. (2.5) 

In addition, we compute /q by examining the fraction of energy that escapes into neutrinos 
in the decay chains of the r lepton. Eqs. (|2.4j ) and (|2.5|) are used to plot the shaded region 
to the right of m x ~ 100 GeV in Figure |2|, namely 



1.2 x 10- 28 -^7 ( GeV ) 5 , for m x < 800 GeV 
s GeV \ m x J 



(2.6) 



2 x 10" 30 ( - — ) , for m v > 800 GeV. 



cm 3 / 100 GeV \ 3 
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For the bremsstrahlung of gamma rays, we take /q = 1. Dark stars can form for m x < 880 
GeV. If /q < 0.86, the bremsstrahlung does not play any role in determining what is 
the largest possible mass of neutralino forming dark stars, and dark stars can form for 
m x < 830 GeV. The annihilation cross section in dark stars can be as low as the lower 
edge of this shaded region, while the correct relic density is obtained through a much larger 
effective annihilation cross section. We notice that in most of the shaded region dark stars 
can still form, except at the higher masses where the shaded region crosses the boundary 
of the area marked 'no dark star.' 

Other coannihilations may arise in the MSSM. For instance, one might have coanni- 
hilations between the neutralino and a sneutrino or a selectron or a smuon. These may 
lead to even smaller (<7i>)ds than the case of stau coannihilations we use as an example, 
and so lead to a situation in which dark stars do not form. The worst case for dark stars 
is coannihilatin with sneutrinos, in that neutrinos are generated in the final state and 
neutrinos escape from the forming protostar without depositing energy. Similarly to the 
neutrinophilic case discussed below, three-body annihilation channels need to be consid- 
ered, in particular internal and final state bremsstrahlung of charged leptons. A simple 
estimate of Z bremsstrahlung in the final state gives us 

„ (trCxX -»• wZ)v) , n cm 3 /100GeV\ 3 . x 

to mx - 3xl ° TGevi^r) ■ (2 - 7 » 



Here we took /g = 1/2 as a representative value. In the process of virtual internal 
bremsstrahung, Z can take away a sizable fraction of the energy and /q may be sizable. 
Eq. ([T?]) is plotted in Figure ^ as the dashed line near the edge of the shaded coannihilation 
region. In terms of dark star formation, this case is similar to coannihilation with the f . 
Dark stars can form up to m x = 1 TeV. For different choice of /q, the cross point is at 
m x = (2/ Q )§ TeV (m x = 800 GeV for f Q = 1/4.) 

We therefore conclude that except in very special cases, namely on top of the Z reso- 
nance or for coannihilations of heavy sleptons or sneutrinos (m x > 800 GeV), dark stars 
can form in the MSSM. 

3. Kaluza-Klein Dark Matter 

If the Standard Model lives in five or six dimensions and the extra dimensions are compact- 
ified at a radius ~ 1 /TeV, the Kaluza-Klein (KK) number can be preserved in a consistent 
way with all the interactions involving an even number of odd KK number particles. In this 
setup, the lightest KK particle (LKP) is stable and can be a good dark matter candidate 
0. In particular, there are two interesting KK candidates for the dark matter: the KK 
photon (more precisely the KK modes of the U(1)y gauge boson) and the KK neutrino. 

The KK photon annihilates to quarks and leptons through the i-channel exchange of 
KK fermions, and its relic abundance is compatible to observation for its mass around 
1 TeV. If the right-handed KK electron, muon, and tau are nearly degenerate (i.e if the 
mass difference is <1%), the KK mass needed for the right relic density can drop to 700 
GeV, since in this case the coannihilation cross section is very small compared to the self 
annihilation cross section, leading to a smaller effective cross section. 

As far as the KK neutrino is concerned, its annihilation cross section to quarks and 
leptons proceeds via t- or s-channel exchange of gauge bosons, while annihilations to gauge 
bosons are mediated by ^-channel KK lepton exchange or s-channel gauge bosons. If one 
flavor of the KK neutrino is considered, the correct relic density is obtained for a mass 
around 1.5 TeV. Including three flavors the effective cross section becomes smaller due to 
coannihilations between different flavors and the mass leading to the correct relic density 
is around 1 TeV. An additional coannihilation process with the KK left-handed electron is 
also possible when the latter has a smaller mass splitting with the KK neutrino, but this 
effect is almost negligible. 

In the case of KK dark matter the s-wave anniilation cross section is always sizable 
both for the KK photon and for the KK neutrino, so that there is little difference between 
{ av )ds an d (<7w)f D . Moreover, the temperature in the dark star is very low compared to 
the freeze out temperature, so coannihilations with other particles give no contributions to 
the effective cross section with the exception of exact degeneracy of the masses. Anyway, 
since for KK dark matter the coannihilation cross section is either smaller than the one 
without coannihilation or the difference between the two is negligible, (crv)ds is expected 
to be always larger or comparable to (av)f . 

Two interesting exceptions to the scenario described above are resonant annihilation 
with level-2 KK particles Q and coannihilation with the KK gluon ||. In principle these 



effects can enhance the effective cross section at the freeze out temperature, so that, if the 
latter is normalized to that of a thermal relic, the cross section in the dark star can be 
suppressed. In particular, the s-channel annihilation at one loop through the exchange of 
the second KK Higgs with mass m h (2) is discussed in ||] . The canonical value for a thermal 
relic (crv)f Q = 3 x 10~ 26 cm 3 s _1 is obtained for a KK photon with mass ttikk — 800 GeV 
if the relation m h (2) = 2?tikk holds up to 5%. This new enhancement changes the cross 
section only by 10 to 20%. Therefore, the largest possible difference between (av)f Q and 
(o"^)ds can be at most 10 to 20%. Unlike the MSSM, there is no p-wave suppression for the 
KK photon annihilation and the cross section at freeze out temperature is large enough due 
to the t-channel exchange of KK fermions. The same is true for KK neutrinos through the 
t-channel exchange of KK Z and KK W bosons. The KK Z and KK W bosons have similar 
masses, and in this case /q <C 1 is not possible. Thus /^(cr^ds/w-KK > 10~ cm 3 s /GeV 
is safely satisfied. Coannihilation with the KK gluon is a last interesting possibility Q. If 
the KK gluon and the KK photon are degenerate with an accuracy much less than 1%, the 
correct relic density is obtained for wikk — 5 TeV. By comparing it to the case without 
coannihilation (tokk = 700 GeV), one can see that coannihilation with the KK gluon can 
enhance the effective cross section by a factor of 50. Even in the worst scenario in which 
coannihilation with the KK gluon is effective at the freeze out temperature but is absent 
in the dark star, /Q^^ds/^KK > 4 x 10 _32 cm 3 s" 1 /GeV (assuming a conservative value 
/q = 1/3) and the dark star can form. 

As a consequence, if (crv){ = 3 x 10 _26 cm 3 s _1 is imposed in order to explain the 
observed relic density, the condition for the dark star formation is always satisfied. We can 
conclude that KK dark matter that explains the observed relic density can always form 
dark stars. 



4. Leptophilic Models 



Leptophilic DM models |1(J have recently become popular in order to explain simultane- 
ously the excess in PAMELA positrons pif| , the excess in Fermi-LAT electrons fl^l , as well 
as the excellent agreement between the observed antiproton spectrum and the correspond- 



ing standard expectation [13]. In order to explain the excesses, the mass of the DM is also 



constrained to be larger than about 100 GeV. This constraint might be more stringent if 
the electron FERMI-Lat data are taken into account, m > 400 GeV. 

In leptophilic models the DM particles generically annihilate exclusively to charged 
leptons, either of only one type (electrons, taus or muons) or democratically to all the 
three families. It is also possible to consider decays to neutrinos, but for simplicity we 
will not consider this case which would simply imply a straightforward generalization (see 
Section 5 for annihilation into neutrinos only). In particular, in this section we discuss the 



case of democratic annihilation to the three lepton families. In this case, using PYTHIA |20 
one gets /q ~ 0.56 almost constant in the range 200 GeV< m < 2 TeV. 

In order to explain the PAMELA and Fermi-LAT excesses, large annihilation cross 
sections 10~ 25 < (<7i>)g a i/cm 3 s -1 < 10 -23 are needed at the velocity of DM particles in 
our Galaxy, v ga i ~ 300 km/s. Assuming (<7-u) ga i = (av)f (s-wave annihilation) these 



values are up to two orders of magnitude larger than the value (crv)f ~ 3 x 1CP 26 cm 3 
s _1 compatible with a standard thermal relic abundance in agreement with observations. 
Clearly, the case of p-wave suppression of (av) ga _i would be even worse. Several mechanisms 
have been devised in order to explain this discrepancy, such as a non-thermal production of 
the DM particles, a non-standard evolution history of the Universe or an enhancement of 
the annihilation cross section at low velocities (Sommerfeld effect). In this sense leptophilic 
models represent another interesting possibility in connection with the formation of a dark 
star: a very large annihilation cross section throughout the history of the Universe and/or 
at the low temperatures where dark stars are formed. 

Assuming (crv} ga i=(av}d s , the previous discussion implies that in a dark star: 

7 x lO-^cm^GeV- 1 < fq^ 1 < 8 * lO^Ws^GeV- 1 , (4.1) 

for 100 GeV < m < 2 TeV. This is shown in Fig. ||, where along with the intervals 
required to explain the PAMELA and Fermi/LAT excesses the present constraints for the 
combination fQ^crvjte/m are summarized as a function of m. In this plot we have assumed 
that (av) does not depend on the temperature, in order to directly compare constraints 
relative to different epochs. 

In particular, the thick and thin solid line closed contours show the range of values 
compatible with the PAMELA positron excess Q and the FERMI-LAT e + + e~ data 



respectively. On the other hand, the two solid open lines represent conservative 2 a C.L. 
upper bounds for fQ(crv)ds/M x obtained from the flux of e + + e~ observed by FERMI [12] 
(thin line) and the e + + e~ flux measured by HESS (l5| (thick line). 

In the same Figure, we also plot with the dotted line the upper bounds on fq{(7v)ds/m 
obtained by comparing the expected gamma-ray flux produced by Inverse Compton (IC) 
scattering of the final state leptons to the diffuse flux of gamma-rays measured by FERMI 



at intermediate Galactic latitudes [17]. 



Finally the long and short dashed line shows the upper bound on fQ{crv)d s /m obtained 
by considering the imprint on the Cosmic Microwave Background Radiation (CMB) from 
the injection of charged leptons from DM annihilations at the recombination epoch Ji~6[ ]. 

It is clear from Fig.||| that the sizable values of the combination fQ{av)^ s /M x are 
compatible to the formation of a dark star (Eq.[Q|) when {av)& s = {<rv) g& \, with the range 



m > 1 TeV disfavored [14] by several constraints. On the other hand, by assuming a 
Sommerfeld enhancement of the annihilation cross section one may have {crv)ds ^ (o"^)gai 
depending on whether the enhancement effect is already saturated at the velocity v > 10 km 
s _1 inside the dark star, possibly implying in this case an even more favorable situation for 
the dark star formation. However, in presence of a non-saturated Sommerfeld enhancement 
at the recombination epoch the CMB constraint could be stronger, since at z ~ 1100 DM 
particles are slower (v ~ 10~ 8 c ~ 10~ 3 km s~ 1 flfj) than inside a dark star. In this case 
leptophilic DM could explain the PAMELA and Fermi/LAT excesses only by assuming a 
boost factor of astrophysical origin such as dumpiness. In any case, barring some specific 
cases such as the presence of resonances in the annihilation cross section associated with 
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Figure 3: Combination fQ{av)d s /m for a DM candidate annihilating democratically to charged 
leptons of the three families as a function of the mass m, assuming that (av) is independent of 
the temperature. The thick solid line contour shows the range of values compatible with the 
PAMELA positron excess jtl); the thin solid contour is the range compatible with the observed 
e + + e~ flux measured by FERMI-LAT j^]; the thin solid open line is the 2a upper bound from 
the observed e + + e~ flux of FERMI; the thick solid open line is the 2a upper bound from the 
e + + e~ HESS measurement ||l5[ ; the dot-dashed line represents the upper bound on (av) from 
CMB Q. The dotted line shows the upper bound on (av) obtained by comparing the expected 
gamma-ray flux produced by Inverse Compton (IC) scattering of the final state leptons to the 
FERMI-LAT measurement of the diffuse gamma ray emission with subtraction of the expected 
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bound states |0| even in this circumstance the bound in Eq. (|T^) would be easily verified 
and a dark star would be formed. 

Thus we can also conclude that leptophilic models also satisfy the condition to form a 
dark star. 



5. Neutrinophilic Models 

As discussed in the previous sections the most popular examples of thermal dark matter 
candidates, namely the neutralino in the MSSM and the KK photon or KK neutrino in 
Kaluza-Klein DM, can easily produce a dark star, provided that the annihilation cross 
section at the temperature of the dark star is similar to that at freeze out and is not 
suppressed by mechanisms such as p-wave annihilation or by the fact that the annihilation 
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cross section is resonant at the freeze out temperature but not inside the dark star. In this 
Section we wish to generalize this statement to the general case of a thermal DM candidate, 
discussing what are the minimal conditions to form a dark star once (crv)ds is normalized 
to the canonical value (crv){ = 3 x 10 -26 cm 3 /s. Moreover, we will also briefly comment on 
the case (crv)^ s < (av)f . 

For this purpose, we need to give a general discussion of the energy fraction fg released 
by DM annihilation into the gas, a quantity that is in general model dependent. Our 
approach to this problem is to consider in this Section the most conservative case of a DM 
candidate annihilating exclusively into neutrinos, i.e. "neutrinophilic" Dark Matter. 

Naively one would expect that the energy fraction of neutrinophilic DM going into 
visible particles vanishes: DM annihilation would generate only neutrinos that would escape 
from the collapsing gas freely. If this were indeed the case, neutrinophilic dark matter 
annihilations would not be able to support a dark star phase. However Z and W bosons 
are expected to be produced from bremsstrahlung radiation of the final state neutrinos, 
so some visible energy, increasing with the mass of the DM particle, is expected to be 
produced by the decay of the Z and / or W. Electroweak bremsstrahlung in the annihilation 
of neutrinophilic DM has already been considered in the context of DM indirect detection 

Let's first assume that neutrinophilic dark matter has an annihilation cross section 
in the dark star equal to the cross section that provides a thermal relic density, namely 
(c^ds = 3 x 10 _26 cm 3 s _1 . At the tree level the branching ratio to vu is 1 if bremsstrahlung 
radiation is neglected. However, when 2m > mw (or mz), on-shell production of W-bosons 
(or Z-bosons) dominates the bremsstrahlung process, which can be viewed as a three body 
decay followed by the subsequent decay of the W or Z gauge bosons. Since the visible 
energy fraction of the W- and Z-boson decays is of order 1, one finds 

/ g 2 E w > g 2 m w 
Wit 2 2m 16tt 2 2m 

Here Ew is the energy of the W boson. In this case, as can be easily checked numerically, 
the condition for forming a dark star is always fulfilled. On the other hand, when 2m < mw, 
off-shell bremsstrahlung occurs, which for 2m <C mw can be treated as a four-body decay 
in the limit of a 4-Fermi interaction. 

The visible energy fraction fg is plotted as a function of the DM mass m in Fig. |2|. 
In this Figure we have used PYTHIA |20| to calculate the subsequent decay of the final 
state particles in the radiative correction, since the value of fg critical for the formation 
of the dark star is near the threshold for the production of an on-shell W-boson, where a 
narrow width approximation or a 4-Fermi interaction are not reliable. From Fig.|| we can 
conclude that a dark star can be formed if the mass of the neutrinophilic dark matter is 
larger than ~ 50 GeV. 

If the DM particle is a scalar or a Majorana fermion, it is possible that the annihilation 
cross section in the dark star is significantly different from the cross section at decoupling. 
A simple example is a scalar neutrinophilic dark matter particle <j) annihilating by t/u- 
channels through the exchange of a heavy fermion. In particular, the s-wave contribution 
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Figure 4: Energy fraction /q released by neutrinophilic dark matter annihilation inside a dark star 
as a function of the DM particle mass m. The black curve represents /q for an ideal neutrinophilic 
DM model where the W and Z bremsstrahlung effects are calculated using PYTHIA p(| . The blue 



dashed line shows the constraint given in Eq. (1.3). Only models above the blue line can form a 
dark star. 

vanishes if the heavy fermion mediating the annihilation is a Dirac fermion and if it interacts 
chirally with the dark matter (j) and the neutrinos, i.e. 

g4>tf L v + My L V R + h.c. + m 2 H 2 , (5.2) 



where M is the mass of the heavy fermion and m is the mass of the dark matter scalar. In 
this case, the annihilation cross section is purely p-wave and is given by 



av 



2 M 4 ™ 2 ((M 2 



m 



2\2 



2m 2 M 2 ) 



(5.3) 



16vr (m 2 + M 2 ) 4 

In this case, since the average velocity of the dark matter in the dark star is ~ 30 km/s, one 
has {<rv)ds = (&v) { (vd s /v{ ) 2 ~ (o"w)fo'yds( m /-^fo) ~ 10~ 32 cm 3 /s. So p-wave annihilating 
neutrinophilic DM cannot make a dark star. 

Notice that in our evaluation of neutrinophilic DM annihilation we have not included 
internal bremsstrahlung of charged particles. Internal bremsstrahlung would increase the 
annihilation cross section of neutrinophilic dark matter, and extend the mass limit for the 
formation of dark stars to values below ~ 50 GeV, but would depend on the specific particle 
physics model. 

Strictly speaking a purely neutrinophilic DM model is not natural, since in most specific 
scenarios other tree-level decay channels contributing to the visible energy are usually 
expected, and a contribution of the latter at the level of ~ 10 -5 level or larger would be 
sufficient to form a dark star. Thus neutrinophilic DM represents a limiting case, allowing 
to show that, as long as {crv)ds — (o'v)fo, any thermal DM candidate heavier than m > 50 
GeV can lead to the formation of a dark star. 
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6. Conclusions 

The first stars to form in the Universe may be powered by the annihilation of weakly in- 
teracting dark matter particles Q. In this paper we explored several popular examples of 
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thermal dark matter models in order to discuss whether they can satisfy the conditions 
for the formation of a dark star: the neutralino in an effective MSSM scenario; leptophilic 
models that might explain recent observations in cosmic rays; the KK-photon and the 
KK-neutrino in UED models; a conservative neutrinophilic model where the dark matter 
particles annihilate exclusively to neutrinos. We find that in general models with thermal 
dark matter lead to the formation of dark stars, with few notable exceptions: heavy neu- 
tralinos in the presence of coannihilations; annihilations that are resonant at dark matter 
freeze-out but not in dark stars; neutrinophilic dark matter lighter than about 50 GeV. In 
particular the discussion of the latter conservative scenario allows us to conclude that a 
thermal DM candidate in standard Cosmology always forms a dark star as long as its mass 
is heavier than ~ 50 GeV and the thermal average of its annihilation cross section is the 
same at the decoupling temperature and during the dark star formation, as for instance in 
the case of a cross section with a non-vanishing s-wave contribution. 

Therefore, we can conclude that the formation of a first generation of stars powered by 
dark matter annihilation is an almost inevitable consequence of thermal dark matter when 
a standard thermal history of the Universe is assumed and if the mechanism of Ref. Q is 
at work. So a dark star is always there whenever there is thermal dark matter. 
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